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Within the frame of kinetic theory a response function is derived for finite Fermi systems which 
includes dissipation in relaxation time approximation and a contribution from additional chaotic pro- 
cesses characterized by the largest Lyapunov exponent. A generalized local density approximation 
is presented including the effect of many particle relaxation and the additional chaotic scattering. 
For small Lyapunov exponents relative to the product of wave vector and Fermi velocity in the sys- 
tem, the largest Lyapunov exponent modifies the response in the same way as the relaxation time. 
Therefore the transport coefficients can be connected with the largest positive Lyapunov exponent 
in the same way as known from the transport theory in relaxation time approximation. 
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The problem of irreversibility is one of the still open 
questions. Two approaches basically can be distinct. One 
approach considers the many particle theory as a suit- 
able starting point to understand the increase of entropy 
as a result of many random collisions leading to irre- 
versible kinetic equations like the Boltzmann equation. 
The other approach considers the theory of determinis- 
tic chaos with the characteristic measure of Lyapunov 
exponent to understand the occurrence of irreversibility. 
While the many particle approach can be easily extend 
to quantum systems the quantum chaos is still a matter 
of debate about the correct term. 

If both approaches describe some facet of irreversibil- 
ity, what we will anticipate in the following, it should 
be possible to give relations between them. While the 
characteristic measure of many body effects is the re- 
laxation time and the transport coefficients, the rele- 
vant measure for chaotic systems is the Lyapunov ex- 
ponent as a measure of phase space spreading of trajec- 
tories. Considerable efforts have been made to connect 
the transport coefficients with the Lyapunov exponent 
In Eld the fact, that the spreading of a small 
phase space volume is given by the sum of Lyapunov 
exponents SV(t) = SV(0) exp Aj)i, is used to give a 
relation between Lyapunov exponents and viscosity. This 
was possible to show with the help of the contact to a heat 
bath in the equation of motion ensuring constant inter- 
nal energy. In the relation between transport coef- 
ficients and Lyapunov exponents was presented in terms 
of Helfand's moments. The interlink was possible to es- 
tablish by reinterpretation of the Helfand's moments as 
stochastic quantities such that the mean variance of the 
time derivatives represents just the transport coefficients. 
In the authors derive a density expansion of largest 
Lyapunov exponent for hard sphere gases from a gener- 
alized Lorentz-Boltzmann equation. This demonstrated 
the intimate relation between transport coefficients and 
dynamical quantities like the Lyapunov exponent. 

Here we like to show that there exists a very simple 



connection between the concept of Lyapunov exponent 
and the dissipation leading to irreversibility for interact- 
ing Fermi systems. It will be shown that if the largest 
positive Lyapunov exponent is smaller than the product 
of Fermi velocity times wavelength in a Fermi system, 
the Lyapunov exponent appears in the same way as the 
relaxation time of the system. Therefore all expressions 
known from kinetic theory, expressing the transport coef- 
ficients in terms of the relaxation time, can be considered 
as an expression in terms of the Lyapunov exponent. 

The concept of response of an interacting many body 
system starts conveniently from the distribution func- 
tion /(p, r, t) satisfying the appropriate kinetic equation, 
which by linearization yields the response to an external 
disturbance. First we discuss the quasiclassical response 
and generalize later to quantum response. The starting 
quasiclassical kinetic equation reads 



9 t /(p,r,i) + -3 r /(p,r,i) 

m 

-a r (Mnd(r, t) + y cxt (r, t))9 p /(p, r, t) 



/o(p,r) - /(p,r,t) 



(1) 



with the self-consistent mean-field potential given as a 
convolution between the two-particle interaction Vq and 
the density Vi n d = / dfVo(r,f)n(r,t), the external dis- 
turbance Vcxt and a typical relaxation time r. The re- 
laxation time approximation serves here as the simplest 
form of collision integral to describe dissipative processes 
by internal collisions of the particles. This leads to a 
natural chaotization and ergodicity of the system. 

Besides this chaotization by mutual collisions we want 
to discuss in the following how additional chaotic pro- 
cesses, e.g. caused by boundary conditions, surfaces etc., 
are influencing the response of the system to external 
perturbation V ext - 

When the equation (|]) is linearized with respect to the 
external perturbation, the selfconsistent potential Vind 
gives a linear density contribution Sn via <5Vi n d = VoSn. 
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Defining the total polarization function as the connection 
between induced density variation and external pertur- 
bation 



with 



5n(x,w) = y<2x'n(x,x',w) SV exi (x', w), 



(2) 



one finds the relation between the polarization function 
including the effect of the selfconsistent potential, II, and 
the polarization without selfconsistent potential, n r , as 

n(x,x') = n r (x,x')+ /dxdxn T (x,x)v (x,x)n(x,x'). 

(3) 

In other words it is sufficient to concentrate on the re- 
sponse function II T to an external potential without self- 
consistent potential T^nd- The selfconsistent response II 
is then given by the solution of the integral equation (||). 
The following derivation of II T is adapted from || . 

Introducing the Lagrange picture by following the tra- 
jectory x(£), p(£) of a particle 

d p 

dt m 

^P = -d x Kxt (4) 
dt 

we linearize the kinetic equation equation (Q) according 
to /(x, p, £) = /o(x, p) + <y/(x, p, £)e _t / T and obtain 

^5/(x(£), P (£), £) = dp/o a x(t) T4 xt . (5) 

This can be integrated to yield 
5/(x,p,£) = 

oo 

-2m /df /dx'A^x' - x(£')) ^4§^Kxt(x', £ + £')• 



<9 P 2 



(6) 



Integrating over p, the density variation caused by vary- 
ing the external potential is obtained as 



/f dp 3 
^ J (27T^j3 9p2 ^ 0( ~ p2 ' X ' 



x f d£'e- l *'^+^V cxt (x', W )-^ ( 5(x'-x(£')) 



(7) 



where g denotes the spin-isospin degeneracy. Comparing 
the expression (T7|) with the definition of the polarization 
function II T in and (|3|) , we are able to identify the po- 
larization of finite systems including the relaxation time 
as 



II T (x,x',w) = n (x,x', uj + -) 

T 



(8) 



f dp 3 

n (x,x',w) = -2mg J j2^Kf dp2 ^°^ 2 ' x ') 
o 

x J d£'e-<'^<S(x'-x(£')). (9) 



Further simplifications are possible if we focus on the 
ground state fo(p 2 ) = 0(p? — P 2 ) of the Fermi system. 
The modulus integration of momentum can be carried 
out and the Kirschnitz-formula P|?[] appears 



n (x,x', = - 



4TT 2 h 3 



<y(x'-x(o)) 



i) 



IU1 J (it I 

— oo 



,//'• / dnp * 



4tt 



<5(x'-x(£')) 



(10) 



This formula represents the ideal free part and a contri- 
bution which arises by the trajectories x(£) averaged over 
the direction at the present time n p Pf = mx(0). In prin- 
ciple, the knowledge of the evolution of all trajectories 
is necessary to evaluate this formula. Molecular dynam- 
ical simulations can perform this task but it requires an 
astronomical amount of memory to store all trajectories. 
Rather, we discuss two approximations which will give us 
more insight into the physical processes behind. First the 
most radical one shows how the local density approxima- 
tion emerges. In the next one we consider the influence 
of chaotic scattering. 

The local density approximation appears from (|1C 
when we perform two simplifications. Introducing 
Wigner coordinates R = (x + x')/2, r = x — x' we have 
to assume 



1. gradient expansion 



P/ (R+-)« P/ (R) + 0(aa) 



(11) 



2. expansion of the trajectories to first order history 

x' - x(£') m -r - £'x + C(£' 2 ) = -r - £'^n p . 

m 

(12) 

With these two assumptions we obtain from ( |l0| ) after 
trivial integrations 

n^ DA (q, R, lu) = J dr e- 4qr n^ DA (r, R, w) 

{oo 
' 

where k = mtv/ (qpf(R)). This can be further integrated 
with the help of 
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dye lky ^-^- = arctan(Im k - iRe k) 1 



2ibx 



l + k 
1 - k 



+ 7r[sgn(l + fc)+sgn(l-fc)] 



Imfe^O 

(14) 



to yield the standard Lindhard result (gTj) in the classical 
limit 



nr(q»p/,w) 



4tt 2 ^ 



1 - k 



+ ikw [sgn (1 + k) + sgn (1 - k)] j (15) 

where & = muj/(qpf). We recognize the ground state re- 
sult for infinite matter except that the Fermi momentum 
p/(R) has to be understood as a local quantity corre- 
sponding to local densities so that we get with rtq) 



nL DA (q,R, w ) = nj ) nf (q,p / (R), 



U + -) 
T 



(16) 



For extensions beyond the local density approximation 
see §§. 

Now we focus on the influence of an additional chaotic 
scattering which will be caused e.g. by a surface bound- 
ary. In order to investigate this effect we add to the 
regular motion (113) a small irregular part Ax 



x' - x(t') « -r - 1' 



Pf 



Ax. 



(17) 



The irregular part of the motion we specify in the direc- 
tion of the current movement lasting a time At and given 
by an exponential increase in phase-space controlled by 
the largest Lyapunov exponent A. Therefore we can as- 
sume [f < 0] 



Ax w ^£A t exp[-A(t' - A t )] + const. 



(18) 



Since we are looking for the largest Lyapunov exponent 
we can take ( |l8| ) at the maximum A t = — 1/A. Further, 
we require, that in the case of vanishing Lyapunov expo- 
nent we should regain the regular motion (|l^). We have 
for (0) therefore 



x' - x(t') n -r - 



1 - exp(-Ai') 



A 



(19) 



With this ansatz one derives from (|T^) instead of Jl3| ) 
the result 



n x (q,R,w) = - 



mgpf(K) 



1 + ik I dy^ I I 



j/X-l 



(20) 



which for A — > resembles exactly (|T3|). The further 
integration could be given in terms of hypergeometric 
functions but this is omitted here. 

With this formula together with (||) and (|]) we have 
derived the main result of a polarization function includ- 
ing the influence of many particle effects and additional 
chaotic process characterized by the Lyapunov exponent 
A. 

For the condition 



A << qvF 



(21) 



with Vf = pf /m the Fermi velocity and q the wave length 
we can use lim (1 + a/x) x = exp(a) under the integral 



of (E0|) and the final integration is performed with the 
result of ( |l6| ) but a complex shift 



n 



\ DA (q, R, w) = (q,p/(R), lo + i(X + -)). 



(22) 



We obtain by this way just the known Matthiessen rule 
which states that the damping mechanisms are additive 
in the damping T = — + A. 

Next we discuss the quantum response function and we 
will see that all discussions outlined above can be straight 
forward applied to the quantum response function. In- 
stead of the quasiclassical kinetic equation (|l]) we start 
now from the quantum kinetic equation fct] 



dtf(p, r, t) + — <9 r /(p, r, t) - \ ( ds ^ P 



U(r-lj\ e^'-P)/(p',r,i) 



/o(p,r) - /(p,r, t) 



(23) 



with U = VSnd + V eK t ■ The gradient expansion in U leads 
to first order the quasiclassical expression ([!]) . We follow 
now exactly the same linearization as above and intro- 
duce the Langrange picture. The trajectories are now 
described instead of (0) by the following set 



d p 
—x = — 

at m 



(24) 



where the arbitrary vector s shows the infinite possibil- 
ities of trajectories by quantum fluctuations. The re- 
sulting polarization function for a finite quantum system 
reads now instead of M) 



n (x,x', u) 

xd. 



9 



ir 2 h 3 J {2nhf 
o 

sin (n s Pf)\ f ^„-it'u, x/ „> 



dp /• ffc sin(|sp) 



di'e^ w 5(x' - x(i') - -). (25) 



Compared with (^) we see that due to quantum fluctua- 
tions an additiona integration s appears. Eq. (Eq) is the 
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quantum generalization of the quasiclassical Kirschnitz 
formula (10) for the response function in finite systems. 

Applying now the same gradient approximation Jl2] ) 
we derive from (|25|) with the help of 



' sin(|sp) g ( sin(±s P/ ) 



r 2 h 3 



0(p 2 f -(p-^ 2 )-®(p 2 f-(P + ^ 2 ) 



(26) 



the quantum Lindhard result 
II (q,R,cj) = 

d P e(p 2 / (R)-(p-M) 2 )_e(p2(R)-( p + M ) 2 ) 



ftpq 



(27) 



in local density approximation. 

The ansatz about additional chaotic processes jl9| ) 
leads then to exactly the same expression (|2^) under the 
condition (21) but with the quantum response J27|) in- 



stead of Uf. 

We like to point out that this result has far reaching 
consequences. With the assumption (21) we have shown 
by this way that the linear response behavior is the same 
if dissipation comes from the relaxation time via collision 
processes in many - particle theories or from the concept 
of chaotic processes characterized by the Lyapunov ex- 
ponent. We can therefore state that for small Lyapunov 
exponent compared to the product of wave length and 
Fermi velocity in a many particle system, the largest Lya- 
punov exponent behaves like the relaxation time in the 
response function. 

Since the transport theory is well worked out to calcu- 
late the transport coefficients in relaxation time approx- 
imation we can express by this way the transport coeffi- 
cients in terms of the Lyapunov exponent alternatively. 
This illustrates the mutual equivalence of the concept of 
Lyapunov exponent and dissipative processes in many- 
particle theories. 

Pavel Lipavsky and Vaclav Spicka are thanked for 
many enlightening discussions and A. Dellafiore for 
bringing the Kirshnitz formula to my attention. 
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